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I.  INTRODUCTION 


Essentially,  the  area  A  (v.e.)iR  under  the  Specific  Heat  (Cp)  versus  Temperature  (T) 
plots  between  temperature  limits,  Experimental  Test  Temperature  (Texp)  and  Reactive 
Temperature  (TR)  is  a  measure  of  how  much  atomic  vibratory  energy  explosives  can 
absorb  before  a  reaction  occurs.  The  reaction  may  be  melting,  phase  change, 
decomposition,  burning,  or  even  detonation.  Thus,  to  a  good  approximation,  it  could  be 
expected  that  if  A(v.e.)  Tr  amount  of  energy  is  suddenly  added  via  impact  shock  loading, 
then  a  reaction  may  occur. 

This  A  (v.e.)  tr  concept,  that  impact  shock  sensitivity  or  shock  induced  reactivity  of 
energetic  materials  could  be  related  to  their  specific  heat  (Cp)  variation  with  temperature, 
was  demonstrated  in  References  [1]  and  [2]  for  RDX,  TETRYL,  PETN,  TNT,  and 
TATB,  which  are  basic  secondary  reactive  compounds. 

In  References  [3]  and  [4],  the  A  (v.e.)jR  ideas  were  demonstrated  for  HMX  and  HNS 
which  are  also  important  basic  secondary  explosive  compounds.  The  impact  shock 
response  of  these  seven  compounds  ranges  from  very  insensitive  to  highly  sensitive. 

Most  of  these  seven  basic  energetic  compounds  have  been  the  main  ingredient  of  useful 
explosive  mixtures. 

One  such  mixture  is  the  plastic  bonded  explosive  designated  as  PBX-9502  that  is 
95  percent  TATB  and  5  percent  KEL-F800  [5,  6].  PBX-9502  has  been  rather  extensively 
tested  via  one-dimensional  shock  loading  at  various  temperatures  [7,  8,  9],  and  its 
thennal  characteristics  have  also  been  experimentally  explored  [10,  11,  12].  Much  of 
this  important  information  has  been  collected,  reported,  or  at  least  mentioned,  in 
Reference  [6].  Some  of  these  data  for  PBX-9502  are  listed  in  Tables  1,  2,  3,  and  4. 

Utilizing  this  information,  A  (v.e.)iR  concept  computations  were  made  for  PBX- 
9502.  The  exploratory  comparative  results  for  this  important  energetic  material  provides 
considerable  affirmation  and  support  for  the  A  (v.e.)jR  theory  of  impact  shock  sensitivity. 

General  details  of  the  exploratory  computations  and  experimental  data  comparisons 
involved  in  a  general  A  (v.c.)tr  assessment  are  contained  in  the  following  sections. 
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II.  ANALYSIS 


This  section  contains  excerpts  from  similar  sections  in  References  [1]  and  [3],  which 
should  be  consulted  for  additional  details. 

For  some  explosives,  a  good  estimate  of  the  critical  particle  velocity,  Upcr  ,  where  a 
reaction  (or  detonation)  occurs  is: 


In  some  circumstances,  a  better  estimate  of  the  critical  particle  velocity  is: 


Where: 

A(v.e.)iR  =  fTR  CpdT  (3) 

JTexp 

=  Thermal  vibratory  energy  per  atom  between  Texp  and  Tr, 

Gram  (Cm/Sec)2. 

Cp  =  Specific  heat  per  atom  as  a  function  of  temperature. 

mAv  =  Average  mass  of  an  atom  in  the  material,  Grams  (Appendix  A). 

Texp  =  Temperature  at  which  experimental  impact  shock  tests  are 

conducted.  This  is  normally  room  temperature  (RT~  300  °K  but 
can  (and  should)  be  done  at  higher  and  lower  temperatures. 

Tr  =  Temperature  at  which  some  thermally  induced  reaction  occurs 
(decomposition,  melting,  phase  change,  detonation,  etc.). 

Upcrr  =  Particle  velocity,  Up,  such  that  the  shock  induced  internal  energy 
(eO  is  equal  to  A(v.e.)TR. 

Upcri  =  Particle  velocity,  Up,  such  that  the  total  shock  induced  energy  (et) 
(kinetic  plus  internal)  is  equal  to  A(v.e.)TR. 

et  =  mAvU2p  =  total  shock  energy  per  average  atom. 


2A(v.e.)iR 


yfl  Upci 
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c,  =  U  =  C|(  =  internal  or  kinetic  energy  of  the  shocked  material  per  average 
atom. 

Some  explosives,  when  heated  to  higher  and  higher  temperatures,  melt  before  they 
explode  (RDX  and  TNT,  for  example)  This  melting  will  require  that  the  heat  of  fusion 
(AHp)  be  absorbed  by  the  material  at  T=Tmelt  conditions  before  the  temperature  will 
increase  [13,14],  Consequently,  if  TEXp  is  less  than  TMelt,  then  the  total  heat  absorbed 
from  T  =  Texp  to  T  =  TEXpl  =  Tr  is: 

A  (v.e.)  T,  =  [  Cpdt  +AHf+  f  CPdT.  (4) 

So,  for  solid  energetic  materials  which  melt  prior  to  explosion,  then  A  (v.e.)iR,  as 
defined  by  Equation  (4)  is  employed  in  Equations  (1)  and  (2)  to  compute  UPcRi  and  UPcR2, 
respectively.  Note  that  melting  is  just  one  example  of  a  phase  transformation  which  may 
require  an  enthalpy  increment  (AHj)  to  be  activated.  For  example,  HMX  can  exist  in 
different  solid  polymorphic  forms.  At  a  certain  temperature,  Tj,  one  form  may  change  to 
another  form  if  the  heat  energy  of  transformation  (A  HT)  is  supplied.  So,  AHT  should  be 
added  to  Equations  (3)  and  (4)  if  TR  is  greater  than  Tt. 

Note  that  A  (v.e.)iR  as  defined  by  Equations  (3)  and  (4)  is  actually  an  enthalpy 
increment  (AH).  However,  it  was  shown  via  numerical  examples  in  Appendix  B  of 
Reference  [1]  that,  under  the  experimental  Cp  acquisition  conditions,  the  pressure  times 

<*Tr  /»Tr 

volume  terms  were  minute  compared  to  the  CP  integral,  CpdT.  Thus,  CpdT  is 

essentially  all  of  the  internal  energy  difference  caused  by  thennal  stimulation  during 
standard  tests  at  atmospheric  pressure  to  detennine  the  specific  heat  characteristics. 

Once  A  (v.e.)TR  and  Upcr  values  are  computed,  the  corresponding  shock  velocity 
(Uscr)  is  ascertained  from  experimental  data  for  Us  as  a  function  of  the  particle  velocity, 
Up.  The  experimental  relationship  is  usually  linear  and  written  empirically  as: 

Us  =  Co  +  S  Up.  (5) 

Table  3  lists  the  constants,  Co  and  S  and  the  data  sources,  for  the  explosives 
considered  in  this  investigation. 

When  Up  =  UPcR  and  Us  =  UsCR  are  determined,  the  shock  pressure  is  computed  from 
the  following  well  known  relation: 

Ps  =  A  Us  Up,  (6) 

3 

Where  p„  =  Material  density  (grams/cm  ). 
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Then  Upcr  ,  Uscr  ,  and  Pscr  may  be  compared  to  experimental  shock-induced 
reaction  threshold  information  to  check  the  validity  of  the  above  A  (v.e.)iR  theory  to 
denote  reactive  conditions  under  impact  shock  stimuli.  The  numerical  computations 
involved  in  a  A  (v.c.)tr  assessment  are  straightforward  and  simple  and  may  be  performed 
with  a  hand-held  calculator. 

It  must  be  emphasized  that  any  possible  effect  of  pressure  on  Cp  is  not  taken  into 
account  in  the  present  analysis.  The  basic  idea  is  that  if  a  quantity  of  thennal  vibration 
energy,  A  (v.e.)TR,  under  quiescent  conditions  is  able  to  create  a  reaction,  then  the  same 
amount  of  energy  added  by  an  impact  shock  (e;  or  et)  should  also  cause  some  type  of 
reaction.  The  shock-induced  reaction  may  not  be  the  same  type  as  the  temperature 
induced  reaction,  but  will  nevertheless  be  a  reaction  of  some  kind.  It  may  be  less  or  more 
severe  than  the  thermally  induced  reaction. 

The  Cp  unit  of  calorie/ (atom  -  °K)  is  employed  in  plots  of  Cp  versus  T  information 

23 

which  are  shown  in  this  report.  This  is  because  the  Boltzman  constant,  kB=0.33  x  10'" 
calories/ (atom  -  °K)  and  the  maximum  Cp  at  high  temperatures  for  many  materials  is 
3ke~  1.0  x  10'"'  calories/ (atom  -  °K).  This  is  a  good  mnemonic  reference  level  for 
comparison  purposes.  It  was  noted  in  Reference  [1]  that  the  average  Cp  per  atom  for 
most  polymers  never  reaches  the  3  kB  level  before  a  reaction  (phase  change,  melting, 
glass-to-rubber  transition,  or  even  detonation)  occurs. 

Actually,  Cp  for  some  atoms  or  combinations  of  atoms,  probably  reaches  the  3  kB 
level  and  causes  a  reaction  at  some  Tr.  But  Cp  for  a  large  number  of  atoms  remains 
much  less  than  2  kg.  Thus,  a  large  amount  of  the  possible  thermal  vibratory  energy  is 
never  activated  and  the  average  Cp  per  atom  remains  relatively  low  [15].  In  many  cases, 
important  temperature-induced  reactions  occur  near  the  average  CP  ~  2  kB  level  at 
moderate  temperatures  (400  to  600  °K). 

It  was  shown  in  References  [1]  and  [3]  that  the  magnitude  of  the  average  Cp  per 
atom  at  a  given  temperature  for  seven  explosives  did  not  differ  very  much  from  each 
other.  The  reaction  temperatures,  TR,  for  these  explosives  did  differ  considerably,  but 
mAv  values  for  these  explosives  were  very  similar.  These  observations,  in  conjunction 
with  Equations  (1)  through  (6),  were  sufficient  to  establish  an  overall  correlation  of 
experimental  “shock  sensitivity”  data  with  reaction  temperatures,  Tr  (primarily  Tmelt 
and  TExpl)- 


This  report  documents  a  similar  correlation  for  PBX-9502  with  Texpl  at  four 
different  TExp  conditions. 
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III.  APPLICATIONS 


Strictly  speaking,  for  a  valid  check  of  the  A(v.e.)xR  concept  as  outlined  in  Section  II, 
the  following  experimental  infonnation  is  necessary  for  the  same  reactive  energetic 
substance. 

Thermal  Properties: 

1 .  Specific  heat  (Cp)  as  a  function  of  temperature  (T). 

2.  Tr  =  Tmelt  (melt),  Tt  (Transformation),  TEXpl  (explode) 

3.  AH  =  Heat  of  fusion  (AHp,  melt)  and  heat  of  transformation  (AHt). 

Impact  Shock  Related  Properties: 

1 .  pa ,  Material  density  (gram/cm  )  at  each  test  condition  temperature,  TEXp. 

2.  Us  (Shock  velocity)  as  a  function  of  particle  or  mass  velocity  (Up)  for  each 
TEXp. 

3.  Threshold  Updt  (or  Psdt  ),  where  below  them  detonation  will  not  occur,  and 
above  them  detonation  will  occur  for  each  TEXp. 


It  is  rather  difficult  to  find  all  of  this  consistent  infonnation  for  each  TEXP.  In  fact,  it 
is  unusual  to  locate  impact  shock  results  for  a  TEXp  other  than  RT  ~  290  -  300°K). 

However,  there  was  one  energetic  material,  PBX-9502,  where  impact  shock-related 
data  as  described  above,  were  available  for  four  TEXP  conditions  (-55  °C,  RT,  75  °C,  and 
252  °C).  Even  so,  some  estimates,  approximations,  and  extrapolations  had  to  be  made  to 
fill  in  certain  data  gaps.  This  is  particularly  true  for  certain  Cp  data  and  the  threshold  Up 
or  Ps  (for  one  TEXP  condition).  Additional  details  of  these  data  extension  procedures  are 
given  in  the  following  subsections. 

A.  PBX-9502  Experimental  Thermal  Data 

Reference  [10]  contains  experimental  specific  heat,  Cp  (Cal/Gram-  °C),  data 
in  the  temperature  range  of  37  °C  (3 10  °K)  to  177  °C  (450  °K).  These  data  had  been 
fit  to  a  straight  line  linear  relation  which  could  easily  be  extended  backwards  to 
TEXp  =  -55  °C.  Table  2  contains  these  data  and  also  lists  their  values  when  converted 
to  Cal/Atom  -  °C  units. 

Figure  1  compares  this  PBX-9502  and  TATB  Cp  [Cal/(Atom  -  °C)  data.  As 
expected,  Cp  for  these  two  materials  does  not  differ  very  much  in  magnitude  because 
PBX-9502  is  95  percent  TATB.  The  TATB  non-linear  (quadratic  fit)  Cp  data  came 
primarily  from  Reference  [5]  with  some  extrapolation  necessary  at  temperatures  higher 
than  300  °C  (573  °K).  This  estimated  extension  of  the  TATB  Cp  to  409  °C  (682  °K)  was 
essential  to  the  A(v.e.)xR  analysis  reported  in  References  [1]  and  [2],  and  the  present 
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report  as  well.  This  was  because  some  estimate  of  the  PBX-9502  Cp  was  needed  up  to  at 
least  396  °C  (699  °K),  which  is  given  in  Reference  [12]  as  an  explosion  temperature 
(Texpl)-  To  acquire  this  estimate,  the  linear  relation  for  the  Cp  of  PBX-9502  was 
extended  from  177  °C  (450  °K)  to  292  °C  (565  °K).  From  that  point,  the  extended  Cp  for 
TATB  was  used  for  PBX-9502  also. 

These  Cp  results  (Cal/(atom  -°C)  for  PBX-9502  are  listed  in  Table  2  and 
displayed  graphically  in  Figure  1 .  The  cross-hatched  areas  shown  below  this  PBX-9502 
Cp  versus  T  plot  (between  TExp  and  Tr  =  TExpl)  are  equal  to  A(v.e.)ER  as  defined  by 
Equation  (3). 

The  authors  of  Reference  [12]  heated  PBX-9502  samples  with  a  High  Energy 
Electron  Beam  (HEEB)  and  provided  the  following  information: 

TExpl  =  Thermal  explosion  temperature 

=  396.0  ±  5.0  °C  =  669.0  ±  5  °K 

=  Temperature  at  which  the  HEEB  heating  rate  equals  the 
chemical  reaction  rate. 

TIT  =  Thennal  ignition  threshold 

=  140.0  ±  3.0  Cal./Gram 

=  Explosive  heat  dose 

=  Amount  of  thermal  energy  the  energetic  material  can 
absorb  before  it  explodes. 

Reference  [12]  also  gave  TEXpl  =  409.0 ±  5.0  °C  and  TIT  =  144.0  ±  5.0 
Cal/Gram  for  TATB.  The  A(v.e.)iR  results  for  TATB  reported  in  Reference  [1]  for  TEXp 
=  RT  and  Tr  =  TEXPL  =  409  °C  agreed  well  with  this  TIT  value. 

The  extrapolations  and  assumptions  involved  in  extending  the  PBX-9502 
experimental  Cp  data  to  lower  and  higher  temperatures  detract  somewhat  from 
credibility,  but  appear  reasonable  compared  to  Cp  data  from  other  explosives.  Also,  as 
detailed  later,  a  certain  A(v.e.)iR  value  compares  very  well  with  the  TIT  value  given 
above  for  PBX-9502.  This  lends  considerable  credence  to  the  PBX-9502 
extrapolations/extensions  as  A(v.e)ER  is  directly  related  to  the  Cp  magnitude  via 
Equation  (3). 


6 


B.  PBX-9502  Experimental  1-D  Impact  Shock  Data 

References  [7]  and  [8]  both  document  shock  sensitivity  (Hugoniot  and  run 
distance  to  detonation)  experimental  data  for  PBX-9502  at  ambient  room  temperature 
(RT  =  TEXp  =12-27  °C).  In  addition,  Reference  [8]  also  contains  shock  sensitivity 
results  for  PBX-9502,  which  was  collected  at  -55  °C,  75  °C,  and  252  °C.  These 
investigators  illustrated  shock  sensitivity  to  test  temperatures  by  comparing  run-distance- 
to-detonation  (XD)  versus  shock  pressure  (Ps)  data  on  double-log  Popalato  (POP)  plots. 
Most  of  these  data  are  also  graphically  displayed  in  Reference  [6], 

From  these  sources,  the  following  information  was  extracted  which  was 
essential  to  the  present  analysis  where  detonation  threshold  particle  velocities  (  Updt  ) 
were  compared  to  theoretical  Upcr  results  defined  in  Section  II: 

a.  p,  (Density,  grams  per  cubic  centimeter)  as  a  function  of  temperature, 
Texp-  These  p„  values  are  listed  in  Tables  3  and  5. 

b.  Us  (Shock  Front  Velocity)  as  a  function  of  Up  (particle  or  mass  velocity) 
generated  in  the  PBX-9502  specimen  for  each  TExp-  The  pertinent  values 
of  Co  and  S,  for  Equation  (5),  are  listed  in  Table  3. 

c.  Threshold  Updt  values,  just  sufficient  to  cause  detonation,  were 
detennined  from  tabular  or  graphical  (POP  plots)  data  for  Ps  (or  Up)  as  a 
function  of  run-distance-to-detonation,  XD,  contained  in  References  [6,  7, 
and  8]. 

The  procedure  to  detennine  the  threshold  UpDT  values  is  illustrated  in  Figures  2 
through  5,  where  Up  versus  XD  is  plotted  on  a  linear  scale  from  the  tabulated/plotted  data 
in  References  [6],  [7],  and  [8]. 

For  Texp  =  RT,  75  °C,  and  252  °C,  tabulated  Up  versus  XD  experimental  data 
asymptotically  approach  a  practical  detonation  threshold  value  of  Updt  as  XD  gets  larger 
and  larger  as  shown  in  Figures  3,  4,  and  5,  respectively. 

The  values  estimated  by  this  procedure  are  0.88,  0.82,  and  0.46  km/sec  for 
Texp  =  RT,  75  °C  and  252  °C,  respectively.  These  estimated  UPDT  threshold  values 
also  appear  reasonable  as  shown  in  Figures  6,  7,  and  8,  where  UP  is  plotted  versus  tD 
(time  to  detonation)  for  TExp  =  RT,  75  °C  and  252  °C,  respectively. 

The  Us  versus  Up  variation  for  PBX-9502  at  RT  conditions  indicates  some  type 
of  impact-induced  reaction  is  occurring  around  Up  =  1.0  ±  0.2  km/sec  where  a  slope 
(dUs/dUp)  change  is  evident  [7,  20],  This  is  consistent  with  the  RT  Updt  estimated  result 
of  0.88  km/sec. 
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The  estimation  of  the  detonation  threshold  Updt  for  Texp  =  -  55  °C  was  not 
quite  so  straightforward  as  that  for  Texp  =  RT,  75  °C,  and  252  °C.  The  reasons  for  this 
are  detailed  as  follows. 

The  shock  sensitivity  data  (Ps,  XD)  from  PBX-9502  experiments  at 
Texp  =  -55  °C  were  apparently  documented  in  Reference  [9].  References  [6]  and  [8] 
show  these  results  in  the  form  of  Xd  versus  Ps  points  on  a  double-log  POP  plot.  For  most 
explosives,  XD  versus  Ps  is  a  linear  (straight  line)  function  on  a  log-log  POP  plot.  This 
seems  to  be  true  for  PBX-9502  at  -55  °C  and  LX- 17  at  ambient  (RT)  conditions,  which 
both  exhibit  the  same  trend  (slope)  and  magnitude. 

However,  the  PBX-9502  data  plot  at  -55  °C  began  at  XD  ~  10mm  with  a 
corresponding  Ps  of  about  140  kbars.  This  is  obviously  somewhat  different  from  the 
threshold  values  for  this  condition.  A  linear  extrapolation  back  to  Xd  =  50.0  mm  gives  a 
Ps  value  of  about  90.0  kbars,  which  should  be  close  to  the  threshold  (go  or  no  go)  value. 
The  corresponding  result  of  UPDT  is  0.95  km/sec,  which  is  found  via  Equations  (5)  and  (6) 
with  the  appropriate  input  conditions. 

Table  4  gives  a  list  of  the  Xd  versus  Ps  points  extracted  from  the  PBX-9502 
data  for  TExp  =  -55  °C,  shown  in  Figure  IV-9  of  Reference  [6].  The  corresponding  values 
of  Up  are  also  listed  in  this  table.  They  are  plotted  versus  Xd  in  Figure  2  of  the  present 
report  where  an  extrapolation  is  made  to  yield  a  threshold  value  of  Updt  equal  to  around 
0.94  or  0.95  km/sec  for  this  TExp  condition. 

All  of  this  extrapolation  to  estimate  this  detonation  threshold  Updt  value  at 
Texp  =  -55  °C  detracts  from  credibility,  but  the  result  seems  reasonable  (with  the  proper 
trend  and  magnitude)  as  shown  in  Figure  9. 

Figure  9  depicts  the  estimated  detonation  threshold  UpDT  values  versus  TExp- 
These  results  are  also  listed  in  Table  5.  It  is  these  Updt  results  (determined  basically 
from  experimental  data)  that  the  computed  Upcri  and  Upcrr  are  compared  with  in  order  to 
check  the  validity  of  the  A(v.e.)TR  concept. 

A(v.e.)iR  magnitudes  (and  thus  Upcri  and  Upcrj  magnitudes)  also  depend 
entirely  on  experimental  information  (Cp  vs.  T  and  Tr).  Therefore,  a  check  on  the 
validity  of  the  A(v.e.)  iR  concept  is  essentially  a  check  on  the  consistency  of  the 
experimental  thermal  and  impact  shock  data. 


C.  PBX-9502  A(v.e.)TR,  UPcR,  and  PSCr  Results 


This  subsection  delineates  the  combination  of  information  from  both  A  and  B 
sub-sections  to  systematically  compute  A(v.e.)  tr,  UpCR,  and  Pscr  results  for  each 
experimental  test  condition  (TEXp)  with  a  plausible  reactive  temperature  (Tr)  value. 

Basically,  A(v.c.)tr  was  computed  via  numerical  integration  of  the  area  under 
the  PBX-9502  CP  versus  T  Data  Plot  (Fig.  1)  between  TEXP  and  Tr.  These  areas  were 
either  trapezoids  or  rectangles  because  the  CP  variation  with  T  was  linear  or  constant  as 
discussed  above. 

A(v.e.)TR  was  evaluated  for  each  of  the  four  temperatures  (-55  °C,  RT,  75  °C 
and  252  °C)  where  experiments  were  conducted.  This  was  done  for  only  one  reactive 
temperature  (TR  =  TEXPL  =  396  °C  =  669  °K)  given  in  Reference  [12]  as  the  temperature 
where  an  explosion  occurred  under  the  rapid  heating  from  a  HEEB. 

PBX-9502  can  explode  at  lower  temperatures,  particularly  if  soaked  for  a  long 
period  of  time  (14  hours)  at  a  constant  temperature  just  above  the  so  called  Arrhenius 
critical  temperature.  This  is  Tcr  =  296  °C  =  542  °K  which  is  described  in  Reference  [11] 
as: 


“That  surface  temperature  at  which  the  internal  energy  generated 
by  chemical  decomposition  is  greater  than  that  which  can  be 
removed  through  the  surface  by  thennal  conduction.” 

If  soaked  at  a  temperature  below  TcR,  apparently  no  explosion  occurs.  TEr  = 
296  °C  =  542  °K  was  not  employed  as  a  reactive  temperature,  Tr,  in  the  present  analysis 
because  the  time  scales  involved  in  these  so  called  Arrhenius  tests  (hours)  and  impact 
shock  loading  (psec)  are  vastly  different.  The  rapid  heating  (energy  input)  during  HEEB 
tests  [12]  is  much  more  compatible  to  impact  shock  loading  energy  deposition  conditions 
than  slow  cook-off  tests. 

Thus,  Tr  =  TEXPL  =  396  °C  =  669  °K  from  the  HEEB  tests  [12]  was  employed 
for  the  A(v.e.)xR  computations. 

Table  6  lists  the  A(v.e.)xR  values  computed  for  each  of  the  four  temperatures 
(Texp)  where  impact  shock  experiments  were  conducted.  Table  6  also  lists  the  associated 
results  for  UPCri  and  UpCr2.  Table  7  lists  the  corresponding  Uscr  and  Pscr  information. 

Upcri  and  UPCR2  are  plotted  versus  Texp  in  Figure  9  where,  for  comparative 
purposes,  the  estimated  detonation  threshold,  UPDT,  results  are  also  shown. 
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IV.  DISCUSSION 


Reference  [8],  in  particular,  illustrated  PBX-9502  impact  shock  sensitivity  to 
different  test  condition  temperatures  (Texp).  This  was  accomplished  via  impact  testing  at 
different  temperatures  and  the  results  were  presented  in  POP  plots  (XD  vs.  Ps)  from 
Figure  5  in  Reference  [8].  There  are  both  amplitude  and  slope  changes  which  illustrate 
impact  shock  sensitivity  to  temperature.  For  instance,  PBX-9502  at  Texp  =  -55  °C  and 
+75  °C  are  slightly  above  and  below,  respectively,  the  POP  plot  for  ambient  RT 
conditions.  So,  there  was  not  a  large  change  in  shock  sensitivity  at  TExp  =  -55  °C  and 
+75  °C  (compared  to  TExp  =  RT),  but  there  was  a  large  change  (increased  shock 
sensitivity)  for  TExp  =  252  °C. 

This  report  attempts  to  provide  a  quantitative  explanation  for  the  temperature- 
related  impact  sensitivity  phenomena  described  above.  These  trends  are  in  agreement 
with  the  detonation  threshold  Updt  data  (Section  III-B)  listed  in  Table  5  and  depicted 
versus  TExp  in  Figure  9.  This  is,  perhaps,  to  be  expected  because  the  Updt  results  were 
derived  from  the  same  experimental  POP  plot  data. 

There  is  not  much  difference  in: 

UpDT  =  0.94  km/sec  TEXp  =  -55  °C 

UpDT:=0.88  km/sec  TExp  =  RT 

UpDT  =  0.82  km/sec  TExp  =  75  °C 

There  is  only  0.06  km/sec  difference  in  Updt  for  ambient  RT  and  Updt  for 
Texp  =  -55  °C  and  Updt  for  TEXp  =  75  °C. 

This  same  trend  shows  up  in  the  computed  results  for  Upcri  and  Upcr2  at 
Texp  =  -55  °C,  RT  and  75  °C.  Fundamentally,  there  is  just  not  that  much  difference 

fTR 

between  CPdT  =  A(v.e.)TR  for  these  conditions,  as  shown  in  Table  7.  Upcr  is 

JTexp 

proportional  to  the  square  root  of  A(v.e.)iR  which  makes  the  differences  in  Upcr  even 
less. 


When  Texp  =  252°,  this  is  much  closer  to  TR  =  396  °C  than  TExp  for  the  other  three 

(•396 

test  conditions.  The  value  of  CpdT  =  A(v.c.)tr  for  this  case  is  less  than  half  of 

J252 

A(v.e.)iR  for  the  other  three  test  conditions  (Table  7).  Consequently,  the  UpCR  values  for 
Texp  =  252  °C  are  considerably  less  than  UpCR  results  for  the  other  three  test  cases  at 
lower  temperatures.  This  is  because  much  more  heat  energy  had  been  added  to  the 
explosive  target  before  it  was  impact  shock  loaded.  This  is  consistent  with  the 
experimental  UpDT  results  for  this  test  case  (Table  5  and  Fig.  9). 
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It  should  be  noted  that  Tr  =  396  °C  for  PBX-9502  is  rather  large  compared  to  most 
explosives  [1  and  3].  For  Tr  values  much  less  than  this,  these  explosives  will  exhibit 
more  variation  in  impact  shock  loading  sensitivity  at  temperatures  closer  to  ambient  RT 
(say  -55  °  or  75  °  C)  than  PBX-9502.  This  follows  from  the  integral  limits  (TEXp  and  TR) 
in  the  definition  of  A(v.c.)tr. 

One  highly  affirmative  check  on  the  A(v.e.)iR  computations  with  TEXp  =  RT  (20  °C) 
and  Tr  =  396  °C  is  as  follows  for  these  PBX-9502  conditions: 

A(v.e.)TR  =  247.8165(10  23)  calories  per  average  atom  (Table  6). 


=  136.0747  cal/gram  because  mAv  =  1.821 18(1  O'23)  grams  per 

average  atom  of  PBX-9502. 

This  compares  rather  well  with: 

TIT  =  140.0  +/-  3.0  cal/gram 

=  thermal  ignition  threshold  from  HEEB  tests  (see  Section  III. A 
and  Reference  [12]). 

This  is  certainly  indicative  of  reasonably  good  values  of  CP  versus  T  information  for 
the  A(v.e.)iR  computation  at  this  condition.  No  adjustments  to  the  Cp  versus  T 
information  were  made  in  order  to  yield  this  good  comparison  between  these  A(v.e.)jR 
and  TIT  values. 
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V.  CONCLUSIONS 


Based  on  the  good  comparative  results  described  in  the  preceding  section  and  the 
previous  work,  documented  in  References  [1]  through  [4],  it  is  concluded  that  the 
A(v.e.)iR  concept  for  explosives: 

1 .  Provides  a  quantitative  connection  between  impact  shock  sensitivity 
(U  Pdt  orPsDi)  and  the  thennal  properties  (Cp  and  Tr). 

2.  Illustrates  the  consistency  between  experimental  impact  shock  results  and  the 
thermal  properties. 

3.  Provides  additional  support  for  the  premise  that  the  specific  heat  per  average 
atom  (at  a  given  temperature  up  to  TR)  has  approximately  the  same  magnitude.  See 
References  [1]  through  [4]  where  this  is  illustrated  for  several  important  secondary 
CHNO  energetic  materials. 

4.  Indicates  that  high  reactive  temperatures  (TR)  are  related  to  phenomenal  impact 
shock  insensitivity  such  as  exhibited  by  TATB  and  PRX-9502.  That  is,  as  TR  becomes 
larger,  higher  impact  shock  pressure  (Ps),  or  actually  larger  particle  velocity  (Up),  is 
required  to  initiate  detonation.  This  is  quantitatively  explained  by  Paragraph  3  above  and 
the  basic  definition  of  A(v.e.)jR  that  is  Equation  (3)  where  the  upper  limit  of  the  integral 
is  Tr.  As  Tr  increases,  so  does  A(v.e.)TR  which  is  the  amount  of  heat  energy  (or  impact 
shock  energy)  that  can  be  absorbed  before  a  reaction  occurs. 
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VI.  RECOMMENDATIONS 


In  order  to  perform  this  comparative  analysis,  the  amount  of  experimental  data 
extrapolation  and  judicious  extension  necessary  to  provide  certain  impact  shock  and 
thennal  information  suggests  that: 

1.  Experimental  CP  versus  T  data  should  be  collected  for  as  wide  a  range  of  T  as 
possible  on  both  sides  of  ambient  (RT)  conditions. 

2.  Plane  impact  shock  detonation  thresholds  (UpDT  of  Psdt)  should  be  delineated  at 
each  temperature  (TExp)  where  these  experiments  are  conducted. 

3.  Along  with  recommendation  No.  2,  Us  versus  Up  data  are  also  needed  at  each 
Texp  in  order  to  compute  UsCr  and  Pscr  from  the  UpCR  results. 

These  are  generalizations  of  the  specific  recommendations  made  in  Reference  [1]. 
They  apply,  of  course,  to  present  explosives,  but  are  particularly  relevant  to  new 
explosives  under  development. 

It  is  realized  that  these  recommendations  may  be  difficult,  or  practically  impossible, 
to  perfonn,  particularly  at  high  temperatures  where  a  violent  reaction  may  occur. 

One  other  recommendation,  or  admonishment,  to  impact  shock  investigators  is: 

Pay  as  much  attention  to  particle  (or  mass) 
velocity  (Up)  magnitudes  as  a  measure  of 
shock  induced  reactions  as  is  currently 
bestowed  on  shock  pressure  (Ps)  levels. 


This  advice  is  solidly  based  on  the  investigations  documented  in  References  [1] 
through  [4]  and  this  report,  plus  additional  studies  concerning  other  shock  loading 
induced  reactions  of  solid  materials  not  cited  herein. 
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Table  1.  Thermal  Properties  of  TBX-95  02 


Property 

Reference 

Value 

Texpl  Thermal  Explosion  Temperature 

[12] 

396.0+/-  5.0  °C 
669.0  +/-  5.0  °K 

TIT  Thermal  Ignition  Threshold 

[12] 

140.0+/- 3.0 
Cal/Gram 

A  HF  Heat  of  Fusion 

[11] 

50.0  Cal/Gram 
at  420.0  °C  (693  °K) 

Cp  Specific  Heat 

[5,  6,  10] 

See  Table  2. 

TCr  (See  IIIC  for  Definition) 

[6] 

[11] 

344.3  °C  to  356.0  °C 

617.3  °K  to  629.0  °K 

269.0  °C 

542.0  °K 
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Table  2.  PBX  -  9502  Specific  Heat  (CP) 


T 

T 

0.00059T(°C) 

Cp 

Cp 

Comments 

°C 

°K 

CAL/(Gram-°C) 

CAL/ (Gram-°C) 

CAL/(Atom-°C) 

— 

-55* 

218* 

-0.03245 

0.21655 

0.39438  *  10-23 

PBX  -  9502  Linear 
Extrapolation 

0 

273 

0.0 

0.24900 

0.45347  *  10-23 

17 

290 

0.01003 

0.25903 

0.47174  *  10-23 

20* 

293* 

0.01180 

0.26080 

0.47496  *  10-23 

27 

300 

0.01593 

0.26493 

0.48249  *  10-23 

37 

310 

0.02183 

0.27083 

0.49323  *  10-23 

PBX  -  9502  Test 
Data  *** 

75* 

348* 

0.04425 

0.29325 

0.53406  *  10-23 

100 

373 

0.05900 

0.30800 

0.56090  *  10-23 

177 

450 

0.10443 

0.35340 

0.64370  *  10-23 

227 

500 

0.13393 

0.38290 

0.69730  *  10-23 

PBX  -  9502  Linear 
Extrapolation 

252* 

525* 

0.14868 

0.39770 

0.72420  *  10-23 

275 

548 

0.16225 

0.41125 

0.74896  *  10-23 

280 

553 

0.16520 

0.41420 

0.75433  *  10-23 

300 

573 

0.75910  *  10-23 

^-End  of  TATB  Data 

396** 

669** 

0.75910  *  10-23 

<-  TATB  Data 
Extrapolation 

*  T  =  TExp  (Experimental  Test  Condition,  [8]) 

**  T  =  Texpl  (PBX  -  9502  Explosion  Temperature  [12]) 

***  PBX  -  9502  CP  =  0.249  +  0.00059T  (°C)  for  37  °C  <  T  <  177  °C,  Ref.  [10] 
PBX -9502  mAv=  1-821 18  *  10~23  Grams/Atom 
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Table  3.  Impact  Shock  Loaded  Material  Data  for  Us  =  Co  +  S  Up  Relation 


Material 

TEXP 

Po 

Co 

S 

Source 

Remarks 

~ 

°C 

Grams 

KM 

~ 

~ 

------- 

°K 

CM3 

SEC 

1 

1 

l 

o 

o 

1.91 

3.31 

1.65 

[8] 

See  Figure  6  in  [8]. 

PBX- 

9502 

218  °K 

t 

t 

1 

1 

1 

20  ±  8  °C 

1.891 

1.857 

3.15 

[7] 

Up  <  0.50  KM/SEC 

1 

293  ±  8  °K 

1 

1 

i 

1 

Non- 

Linear 

See 

Remarks 
&  Fig.  4 
in  [71 

Us  =  1.392201  + 
5.153578  Up 

1 

RT 

[7] 

-2.421567  Up2 
+0.561615  Up3 

1 

t 

For  0.5  <  Up  <  1.2 
KM/SEC 

1 

t 

i 

t 

2.938 

1.77 

[7] 

1 .2  <  Up  <2.3 
KM/SEC 

1 

1 

1 

RT 

1.891 

1.90 

3.00 

[20] 

Up  <  .082  km/sec 

1 

1 

1.891 

2.90 

1.78 

[20] 

Up  >  0.82  km/sec 

1 

75  ±  2  °C 

1.857 

2.60 

1.91 

[8] 

See  Figure  6  in  [8]. 

1 

348  ±  2  °K 

1 

t 

1 

1 

1 

252  ±  2  °C 

1.70 

1.33 

3.08 

[8] 

See  Figure  6  in  [8]. 

1 

525  ±  2  °K 

1 

t 
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Table  4.  PBX-9502  Up  Versus  Xd  Information  From  Popolato  Plots  in 
References  [6]  and  [8]  for  TExp  =  -55  °C  =  218  °K 


XD 

[6,8] 

Ps 

[6,8] 

Comment 

po 

Up 

Us 

Ps 

MM 

KBAR 

GRAM/CC 

KM/SEC 

KM/SEC 

KBAR 

1.4 

250.0 

Note  1 

1.91 

2.00 

6.6100 

252.50 

mm 

B a 

4/ 

1.88 

6.4120 

230.24 

y 

4/ 

1.74 

6.1810 

205.42 

sa 

4/ 

1.47 

5.7355 

161.04 

4/ 

1.33 

5.5045 

139.83 

4/ 

1.19 

5.2735 

119.86 

■ 

4/ 

4/ 

1.08 

5.0920 

105.04 

4/ 

4/ 

1.00 

4.9600 

94.74 

4/ 

4/ 

0.95 

4.8775 

88.50 

■ 

4/ 

4/ 

0.94 

4.8610 

87.27 

Note  1.  From  Figure  IV-9  in  [6]. 
Note  2.  Linear  Extrapolation 
Note  3.  Us  =3.31  +  1.65UP 

Note  4.  P§  =  Po  Us  Up 
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Table  5.  PBX-9502  Estimated  Impact  Shock  Induced  Threshold  Infonnation 

From  Up,  Vs,  XD  Data 


Texp 

Texp 

XD 

UpDT 

USdt 

Psdt 

Comments 

□ 

°C 

°K 

MM 

KM/SEC 

KM/SEC 

KBA 

R 

GR/CC 

-55 

218 

19 

4.8775 

4.8610 

88.50 

87.27 

1.91 

See  Figure 

2. 

12- 

28 

(RT) 

285- 

301 

(RT) 

-40.0 

39.9 

1 

4.4348 

4.46 

EXP. 

73.92 

74.80 

EXP. 

1.891 

See 

Figures  3 
&  6. 

75 

348 

-40.0 

0.82 

4.1662 

63.44 

1.857 

See 

Figures  4 
&  7. 

252 

525 

-40.0 

0.46 

2.7468 

21.48 

1.70 

See 

Figures  5 
&  8. 
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Table  6.  Computation  of  UpCRi  and  UpCR2 


ITEM 

mAV 

Texp 

Tr 

Remarks 

/•Tr 

f  CpdT 

JTexp 

(•Tr 

f  CpdT 

JTexp 

(•Tr 

f  CpdT 

JTexp 

AHt 

A(v.e.)TR 

A(v.e.)TR 

UpcRl 

UpCR2 

mAV 

~ 

Grams 

*10-23 

°C 

°K 

°c 

°K 

~ 

Cals 

*10'23 

Joules 

*io-23 

G(Cm/Sec)2 

*10'13 

G(Cm/Sec)2 

*10'13 

G(Cm/Sec)2 

*10'13 

(Cm/Sec)2 

*jg+10 

Km/Sec 

Km/Sec 

PBX- 

9502 

1.82118 

-55 

218 

396 

669 

(Texpl) 

No  AH 

280.4168 

1,173.264 

1.173264 

0.00 

1.173264 

0.64423 

0.8026 

1.1351 

PBX- 

9502 

1.82118 

20 

293 

RT 

396 

669 

(Texpl) 

No  AH 

247.8165 

1,036.864 

1.036864 

0.00 

1.036864 

0.56933 

0.7545 

1.0671 

PBX- 

9502 

1.82118 

75 

348 

396 

669 

(Texpl) 

No  AH 

220.0683 

920.7658 

0.9207658 

0.00 

0.9207658 

0.50559 

0.7110 

1.0056 

PBX- 

9502 

1.82118 

252 

525 

396 

669 

(Texpl) 

No  AH 

108.7080 

454.8343 

0.4548343 

0.00 

0.4548343 

0.24975 

0.4997 

0.7067 
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Table  7.  Computation  of  UsCR  and  PsCR 


ITEM 

Texp 
°C 
°  K 

Po 

Grams/Cm3 

Co 

Km/Sec 

S 

Remarks 

UpcRi 

Km/Sec 

UsCRi 

Km/Sec 

PsCR, 

Kbars 

tr 

°C 

°  K 

UpcR2 

Km/Sec 

UScm 

Km/Sec 

PScM 

Kbars 

PBX-9502 

-55 

218 

1.910 

3.31 

1.65 

No  AH 

0.8026 

4.6343 

71.042 

396 

669 

(Texpl) 

1.1351 

5.1351 

112.368 

PBX-9502 

20 

293 

(RT) 

1.891 

Non  -  Linear 
See  Table  3. 

No  AH 

0.7545 

4.1433 

59.115 

396 

669 

(Texpl) 

1.0671 

4.8166 

97.193 

PBX-9502 

75 

348 

1.857 

2.60 

1.91 

No  AH 

0.7110 

3.9580 

52.259 

396 

669 

(Texpl) 

1.0056 

4.5207 

84.419 

PBX-9502 

252 

525 

1.700 

1.33 

3.08 

No  AH 

0.4997 

2.8691 

24.373 

396 

669 

(Texpl) 

0.7067 

3.5066 

42.128 

20 


1 


kB  =  0.33  (10-23)  CAL/(ATOM-  °K) 


1  EXP 

-55°  C 


1  EXP 

20°C 


1  EXP 

75°  C 


1  EXP 

252°  C 


"*"r  "  "*"expl 
396°  C 


218°  K 


293°  K 


348°  K 


525°  K 


669°  K 


Figure  1.  Specific  Heat  per  Average  Atom  for  PBX-9502,  TATB,  and  RDX 
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(Km  /  sec) 


XD  (Millimeters) 

Figure  2.  U P  Versus  XDfor  PBX-9502  at  TEXP  =  - 
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55  °C  =  218  °K 


Xd  (Millimeters) 

Figure  3.  U  Versus  Xp  for  PBX-9502  at  7A  v/,  =  Ambient  Room  Temperature  (R  T) 
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Figure  4.  UP  Versus  XDfor  PBX-9502  at  TEXP  =  75  °C  =  348  °K 
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Figure  5.  UP  Versus  XDfor  PBX-9502  at  TEXP  =  252  °C  =  525  °K 
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Figure  6.  UP  Versus  tDfor  PBX-95 02  at  Texp  =  RT  Conditions 
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tD  (Micro-Seconds) 

Figure  8.  UP  Versus  tDfor  PBX-9502  at  TEXP  =  252  °C  =  525  "K 
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APPENDIX  A 


Computation  of  mAv  and  diav  for  PBX-9502 

The  solid  materials  considered  in  this  study  were  chemical  mixtures.  For  these  mixtures, 
the  weighted  average  mass,  mAv,  of  a  single  atom  in  the  material  was  desired. 

First,  it  was  necessary  to  compute  the  mass  of  a  single  atom  for  each  of  the  elements 
contained  in  the  solid.  Each  solid  was  composed  of  one  or  more  of  the  following 
elements: 

Carbon,  C;  Hydrogen,  H;  Nitrogen,  N;  Oxygen,  O;  Chlorine,  Cl;  Flourine,  F. 

The  mass  of  a  single  atom  of  these  elements  is: 

MW  gram/(gram-  mole)  _  grams 
Nav  atoms/(gram-mole)  atom 


Where: 


MW= 


grams 
gram  -  mole 


n  f/XtYl  Q 

Nav  =  Avagadros  Number  =  6.02252  x  10+2~'  - 

gram  -  mole 

Table  A-l  lists  MW  and  m  for  each  of  the  elements  in  the  above  list.  Values  of  Nav  and 
MW  are  from  various  chemistry  text  books  and  handbooks. 

To  compute  the  average  weight  (mAv)  of  an  atom  in  the  material,  the  chemical  formula  or 
proportional  chemical  composition  is  required.  Of  course,  the  weight  (m)  of  each 
elemental  atom  must  be  known,  since  niAv  is  just  a  weighted  average  of  the  elemental 
atoms  in  the  material.  The  procedure  is  valid  for  mixtures  of  compounds  as  well  as 
compounds.  See  References  [1-3]  for  examples. 

When  mAv  is  computed,  then  the  average  space  between  the  atoms  (diAv)  is  given  by  the 
following  relation: 

diav  =  (  —  f=cm  [A-2] 

P 

Computations  of  mAv  and  diAv  for  PBX-9502  are  included  in  this  appendix. 


A-l 


Table  A-l.  Mass  of  a  Single  Atom  for  Selected  Elements 


Element 

MW 
Grams 
Gram  -  Mole 

Nav 

Atoms 
Gram  -  Mole 

m 

Grams 

Atom 

Carbon,  C 

12.011 

6.02252(1 023) 

1.9943(1  O'23) 

Hydrogen,  H 

1.008 

6.02252(1 023) 

0.1 674(1  O'23) 

Nitrogen,  N 

14.008 

6.02252(1023) 

2.3259(  1 0"23) 

Oxygen,  0 

16.00 

6.02252(1 023) 

2.6567(10'23) 

Chlorine,  Cl 

35.45 

6.02252(1 023) 

5.8874(10'23) 

Flourine,  F 

18.9984 

6.02252(1 023) 

3. 1546(1  O'23) 

A-2 


Table  A-2.  PBX-9502  itiav  and  diAv 


Chemical  Composition:  C2.30H2.23N2.21  O2.21  Clp^gF 


0.13 


[5] 


C2.30: 

2.30  x  1.9943(10'23) 

— 

4.58689(10  "3) 

GRAMS 

H2.23 

2.23  x  0.1 674(1  O'23) 

= 

0.37330(10'23) 

H2.21 

2.21  x  2.3259(10'23) 

= 

5. 14024(1  O'23) 

02.21 

2.21  x  2.6567(10'23) 

= 

5.87131(1 0"23) 

Clo.038 

0.038  x  5.8874(10'23) 

= 

0.22372(1  O'23) 

Fo.13 

0.13  x  3.1546(10'23) 

= 

0.4 10048(1  O'23) 

9.118  ATOMS 

16.60560(10'23)  GRAMS 

GRAMS 

mAV 

= 

= 

1.821 18(1  O'23)  GRAMS/ATOM 

For  p0 


1.891  GRAMS/CC 


di 


AV 


m,v 


18.21 18(10'24) 
1.891 


9.6308(1  O'24)  CM3 


di 


AV 


2.12759(10'8)  CM 


2.1276  A 


A-3/(A-4  Blank) 
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Appendix  B 

Specific  Heat  (Cp)  Conversion  Relations 

Experimental  Specific  heat  data  are  usually  given  in  calories  (cal)  or  Joules  (J)  per  molecular 
weight  (MW)  per  °  centigrade  (°C)  or  °  Kelvin  (°K).  That  is: 

Cp  =  Cp  (cal.  per  MW  per  °  C  or  °K) 

or 

Cp  =  Cp  (J  per  MW  per  °C  or  °K). 

Occasionally,  Cp  is  presented  as  cal.  or  J  per  gram  °C  or  °K. 

Cp  =  Cp  (cal  per  gram  per  °C  or  °K) 
or 

Cp  =  Cp  (J  per  gram  per  °C  or  °K) 

Note  that: 

Cp  (cal.  or  J  per  gram  per  °C  of  °K)  =  Cp  (cal,  or  J  per  MW  per  °C  or  °K)/MW. 

Since  interest  is  in  Cp  (cal  or  J  per  atom  per  °C  or  °K)  then,  Cp  (cal.  or  J  per  atom  per  °C  or  °K 
mAv  Cp  (cal.  or  J  per  gram  per  °C  or  °K). 

See  Appendix  A  for  procedures  to  determine  mAv-  Note  that  1°  C  =  1°  K,  and  that  the 
conversion  factor  between  Joules  and  calories  is  4. 184  so  that: 

Cp  (J  per  atom  per  °K)  =  4. 184  *  Cp  (cal  per  atom  per  °K). 

Also,  since  1  Joule  =  107  ergs  =  107  grams  (cm/sec)  2 

then, 

Cp  grams  (cm2/sec2)  per  atom  per  °K  =  107.  (J  per  atom  per  °K). 

The  above  relations  were  employed  in  Reference  2,  and  the  present  report.  However,  Cp  can 
also  be  given  in  terms  of  a  velocity  squared  per  °K  as  follows: 

2  2  7 

Cp(cm  /sec"  per  °K)  =10  x  Cp  (J  per  gram  per  °K).  In  these  units  for  Cp,  A(v.e.)iR  = 
f  CpdT  =  (cm/sec)2,  so  Upcri  =  JA(v.c.)tr  =  and  Upcrj  =  J 2A(v.e)iR  =  -\/2  Upcri  . 

JTexp  v  v 


B-l/(B-2  Blank) 
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